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ABSTRACT 

We present statistical analysis of diffuse Galactic synchrotron emission and polarized thermal emission from 
dust. Both Galactic synchrotron emission and polarized thermal emission from dust reflect statistics of magnetic 
field fluctuations and, therefore. Galactic turbulence. We mainly focus on the relation between observed angular 
spectra and underlying turbulence statistics. Our major findings are as follows. First, we find that magnetohy- 
drodynamic (MHD) turbulence in the Galaxy can indeed explain diffuse synchrotron emission from high galactic 
latitude. Our model calculation suggests that either a one-component extended halo model or a two-component 
model, an extended halo component (scale height > Ikpc) plus a local component, can explain the observed an- 
gular spectrum of the synchrotron emission. However, discrete sources seem to dominate the spectrum for regions 
near the Galactic plane. Second, we study how star-light polarization is related with polarized emission from 
thermal dust. We also discuss the expected angular spectrum of polarized emission from thermal dust. Our model 
calculations suggest that C/ cx Z""/^ for / > 1000 and a shallower spectrum for / < 1000. 

Subject headings: MHD — turbulence — ISM:general — cosmic microwave background — Galaxy: structure 



1. INTRODUCTION 

Diffuse Galactic synchrotron emission is an important fore- 
ground source of the cosmic microwave background (CMB) 
signals. Therefore proper understanding of Galactic syn- 
chrotron emission is essential for CMB studies. In this pa- 
per, we present statistical analysis of a synchrotron foreground 
emission map. In particular, we provide physical interpretation 
of the angular spectrum of the synchrotron intensity. In addi- 
tion, we briefly discuss the properties of a model dust emission 
map, another important ingredient of the CMB foreground. As 
measurements of the polarized CMB signals become possible, 
more accurate removal of Galactic polarized foreground emis- 
sion is required. Synchrotron and dust emissions are known 
to be the most important sources of polarized foreground ra- 
diation. Therefore, measurements of angular power spectra of 
such foregrounds are of great interest. In this paper, we also 
provide estimation of angular spectrum of polarized emission 
by foreground dust. 

The angular spectrum of the synchrotron emission delivers 
valuable information on the structure of the Galaxy. The ob- 
served spectra of synchrotron emission and synchrotron po- 
larization (see papers in de Oliveira-Costa & Tegmark 1999) 
reveal a range of power-laws. Since the Galactic synchrotron 
emissivity is roughly proportional to the magnetic energy den- 
sity 0, angular spectrum of synchrotron emission reflects statis- 
tics of magnetic field fluctuations in the Galaxy (see 33.1l for 
further discussions). 

The interstellar medium (ISM) is turbulent and Kolmogorov- 
type spectra were reported on the scales from several AUs to 
several kpc (see Armstrong, Rickett, & Spangler 1995; Lazar- 
ian & Pogosyan 2000; Lazarian 2009). It is believed that mag- 
netic field lines are twisted and bend by turbulent motions in 
the Galaxy. Therefore it is natural to think of the turbulence as 



the origin of the magnetic field fluctuations and thus the diffuse 
synchrotron foreground radiation. Indeed, several earlier stud- 
ies addressed this issue. Tegmark et al. (2000) suggested that 
the spectra may be relevant to Kolmogorov turbulence. Chep- 
urnov (1999) and Cho & Lazarian (2002; hereinafter CL02) 
used different approaches, but both showed that the angular 
spectrum of synchrotron emission reveals Kolmogorov spec- 
trum (C/ oc /""/^) for large values of multipole /. However, they 
noted that the spectrum can be shallower than the Kolmogorov 
one for intermediate values of multipole I, due to density strat- 
ification in the halo (Chepurnov 1999) or the Galactic disk ge- 
ometry (CL02). In this paper, we further elucidate the relation 
between 3 -dimensional turbulence spectrum and observed an- 
gular spectrum. We also investigate how structure of the Galac- 
tic halo affects observed angular spectrum. 

Thermal emission from dust is also an important source 
of foreground emission. Dust emission is the most pro- 
nounced emission in far infrared (FIR) wavelengths. Schlegel, 
Finkbeiner, & Davis (1998) combined 100/im maps of IRAS 
(Infrared Astronomy Satellite) and DIRBE (Difuse Infrared 
Background Experiment on board the COBE satellite) and re- 
moved the zodiacal foreground and point sources to construct 
a full-sky map. Finkbeiner, Davis, & Schlegel (1999) extrapo- 
lated the 100/im emission map and 100/240/im flux ratio maps 
to sub-millimeter and microwave wavelengths. In this paper, 
we present statistical analysis of a model dust emission map 
above. 

Thermal radiation from dust becomes polarized when dust 
grains are aligned. There is ample evidence that dust grains are 
aligned with respect to magnetic field. Therefore, polarization 
by dust is also related to the magnetic field fluctuations. Then 
how are they related? Polarized dust emission is difficult to 
observe directly. Therefore in this paper we first study the re- 



' Dept. of Astronomy and Space Science, Chungnam National Univ., Daejeon, Korea; cho@canopus.cnu.ac.kr 
^ Dept. of Astronomy, Univ. of Wisconsin, Madison, WI53706, USA; cho@..., lazarian@astro.wisc.edu 

^ In fact, correct scaling is e(r) oc n(e)\B± where e(r) is the synchrotron emissivity, n{e) is the cosmic-ray electron number density, and B± is the strength of the 
magnetic field component perpendicular to the line of sight. In our Galaxy, we have 7 ~ 2. Therefore, when is more or less uniform, we have e(r) oc \Bj_ \-. 

1 



2 



Cho & Lazarian 



lation between star-light polarization and turbulence statistics. 
Then we describe how we infer the relation between star-light 
polarization and polarized thermal emission from dust. 

In this paper, we show how MHD turbulence is related with 
the observed angular spectra of Galactic synchrotron emission 
and starlight polarization. For this purpose, we use an analyti- 
cal insight obtained in Lazarian (1992, 1995ab) and numerical 
results obtained in CL02. We also discuss how we can estimate 
polarized microwave dust emission using star-light polarimetry. 
This problem is of great importance in view of recent interest 
to the foregrounds to the CMB polarization. In §2, we briefly 
describe the data sets we use in this paper, the 408MHz Haslam 
map and a model dust emission map at 94GHz. In §3, we re- 
view a simple model of the angular spectrum of synchrotron 
emission arising from homogeneous MHD turbulence. In §4, 
we present statistical analysis of the Haslam map, which is 
dominated by diffuse Galactic synchrotron emission. In §5, we 
investigate the dust emission map. In §6, we consider polarized 
emission from thermal dust. In §7, we discuss how to utilize 
our results to remove Galactic foregrounds. In §8, we calcu- 
late high-order structure functions of the synchrotron and the 
dust maps and we compare the results with those of turbulence. 
Finally in §9, we give summary. 

2. DATA SETS 

We use the 408MHz Haslam all-sky map (Haslam et 
al. 1982) and a model dust emission map that are available on 
the NASA's LAMBDA website 0. Both maps were reprocessed 
for HEALPix (Gorski et al. 2005) with nside=512. 

The original Haslam data were produced by merging several 
different data-sets. "The original data were processed in both 
the Fourier and spatial domains to mitigate baseline striping and 
strong point sources" (see the website for details). The angular 
resolution of the original Haslam map is ^ 1°. Galactic dif- 
fuse synchrotron emission is the dominant source of emission 
at 408MHz. 

The 94 GHz dust map is based on fits to data from earlier 
100 micron and 100/240 micron maps (Schlegel et al. 1998) 
and extension to COBE/FIRAS frequencies and identical to the 
two-component model 8 of Finkbeiner et al. (1999). 

3. SPECTRUM AND STRUCTURE FUNCTION OF DIFFUSE 
SYNCHROTRON EMISSION: A MODEL FOR 
HOMOGENEOUS TURBULENCE REVISITED 

Suppose that 3-dimensional MHD turbulence has a 3D spa- 
tial power spectrum of the form £30 cx k~"\ where k is the 
wavenumber. Note that in Kolmogorov turbulence m= 11/3. 
Then what will be the 2-dimensional angular spectrum, C/, of 
the observed synchrotron total intensity? We mostly follow dis- 
cussions in Cho & Lazarian (CL02). Although we focus on 
synchrotron emission here, the discussion in this section can 
be applicable to any kind of emission from an optically thin 
medium. 

3.1. MHD turbulence and synchrotron emission 

For synchrotron radiation, emissivity at a point r is given by 
e(r) oc n{e)\B±X' , where n(e) is the electron number density, B±^ 
is the component of magnetic field perpendicular to the line of 
sight. The index 7 is approximately 2 for radio synchrotron 
frequencies (see Smoot 1999). If electrons are uniformly dis- 
tributed over the scales of magnetic field inhomogeneities, the 

,http://Iambda.gsfc.nasa.gov/| 



spectrum of synchrotron intensity reflects the statistics of mag- 
netic field. For small amplitude perturbations {6b/B<^\; this is 
true for scales several times smaller than the outer scale of tur- 
bulence if we interpret B as local mean magnetic field strength 
and 5b as random fluctuating field in the local region), if 5b has 
a power-law behavior, the synchrotron emissivity will have the 
same power-law behavior (see Getmantsev 1959; Lazarian & 
Shutenkov 1990; Chepurnov 1999). Therefore, we expect that 
the angular spectrum of synchrotron intensity also reflects the 
spectrum of 3-dimensional MHD turbulence. 

When an observer is located inside a turbulent medium, the 
angular correlation function, hence the power spectrum, shows 
two asymptotic behaviors. When the angle is larger than a crit- 
ical angle, we can show that the angular correlation shows a 
universal 6'"' scaling. On the contrary, when the angular sepa- 
ration is smaller than the critical angle, the angular correlation 
reflects statistics of turbulence. In this small angle limit, we can 
show that the angular power spectrum is very similar to that of 
turbulence. The critical angle is determined by the geometry. 
Let the outer scale of turbulence be L and the distance to the 
farthest eddies be d„,a^. Then the critical angle is 

e^L/d„,a.. (1) 

3.2. Small-angle limit in homogeneous turbulence 

When the angle between the lines of sight is small (i.e. 6 < 
L/d„,ax), the angular spectrum C/ has the same slope as the 
3-dimensional energy spectrum of turbulence. Lazarian & 
Shutenkov (1990) showed that if the 3D spatial spectrum of 
a variable follows a power-law, E^oik) oc k~", then the 2- 
dimensional spectrum of the variable projected on the sky also 
follows the same power-law, 

Ci cx r"' (2) 
in the small 9 limit. For Kolmogorov turbulence (£3/5 oc A;""/^), 
we expect 

c,«r"/3, ife<L/d„,a,. (3) 

Note that I ^tt/O. 

In some cases, when we have data with incomplete sky cov- 
erage, we need to infer C/ from the observation of the angular 
correlation function 

K{e)=<T(ti)T{i2)>, (4) 
where Ci and 62 are unit vectors along the lines of sight, 6 is the 
angle between Ci and 62, and the angle brackets denote aver- 
age taken over the observed region. As we discuss in Appendix 
A, when the underlying 3D turbulence spectrum is oc A;"™ (e.g. 
m= 11/3 for Kolmogorov turbulence), the angular correlation 
function K{9) is given by 

KiO) oc < > -const 6l'"-^ if 6* < L/£/„„^. (5) 
It is sometimes inconvenient to use the angular correlation func- 
tion in practice to study turbulence statistics because of the con- 
stant <T^ >. 

A better quantity in small-angle limit would be the second- 
order angular structure function: 

D2(0) = < |r(ei)-r(e2)P> (6) 

= 2<T^> -2K(e). (7) 
Thus, in homogeneous turbulence with 3D spatial spectrum of 
E(k) oc k~'", we have 

D2i9) oc 9'"-^. (8) 
When we measure the slope of the angular structure function, 
we can infer the slope of the 3D spatial power spectrum of tur- 
bulence. 
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Fig. 1. — Two limits in homogeneous turbulence. Upper plot: Small d limit (9 < L/dmax)- The fluctuations along the entire length of the lines of sight are 
correlated. Lower plot: Large 9 limit (9 > L/dmax}- Only points close to the observer are correlated. Note the definition of r and tj). From CL02. 



3.3. Large-angle limit in homogeneous turbulence 

In this Hmit, the angular correlation function is more useful 
than the structure function. Following Lazarian & Shutenkov 
(1990), we can show that the angular correlation function for 

9 > L/dmax follows 



K(e) = J J dridr2JC(\ri-r2\), 



3 fTT-e/i 
dr rlC(r) / dip oc 

J9/2 



const 



&m6 



(9) 



where we change variables: {r\,r2) ir,'^), which is clear 
from Fig.[T] We accounted for the Jacobian of which is r/sin^. 
We can understand \/9 behavior qualitatively as follows. When 
the angle is large, points along of the lines-of-sight near the 
observer are still correlated. These points extend from the ob- 
server over the distance oc 1/ sin(6'/2). 

If we assume L/d,„ax < 6* <C 1, we can get the angular power 
spectrum C/ using Fourier transform: 



Q 



dOxdOy 



dO 9Jo{W)Ki0) oc r 



(10) 



where 9 = {9]. + 9^)^l^, Jq is the Bessel function, and we use 

K{9)o^9-\ 

3.4. Expectations for homogeneous turbulence 

In summary, for homogeneous Kolmogorov turbulence, we 
expect from equations (O and ( fTOb that 



(11) 



which means that the power index a of C/ i^-l <a< —11/3. 
For small-angle limit, we expect the following scaling for the 
second-order angular structure function: 

D2{9)o,9'l^ if 9 <L/d„,ax. (12) 

The critical angle 9cr ^ L/d,„ax depends on the size of the large 
turbulent eddies and on the length of the line of sight. If we 

^ Note that point sources would result in a 0. 



assume that turbulence is homogeneous along the lines of sight 
and has L ~ 100 ;?c corresponding to a typical size of the super- 
nova remnant, and that dmaK ^ 1 kpc for synchrotron halo (see 
Smoot 1999), we get 9cr 6°. 

4. PROPERTIES OF DIFFUSE GALACTIC SYNCHROTRON 
EMISSION 

In this section, we analyze the Haslam 408MHz all-sky map, 
which is dominated by Galactic diffuse synchrotron emission. 
Our main goal is to explain the observed synchrotron angular 
spectrum using simple turbulence models. 

4. 1 . General properties of the 408MHz Haslam map 

Fig. |2] shows average intensity as a function of galactic lat- 
itude b. We only use data for the Galactic Northern sky. The 
synchrotron emission is roughly constant for high galactic lati- 
tudes (b > 30°), as noted by earlier studies (e.g. de Oliveira- 
Costa et al. 2003). The usual average intensity (solid line) 
and the average taken near the Galactic center show a sharp 
rise towards the Galactic plane. But, the average taken over 
the Galactic anti-center does not show such a sharp increase to- 
wards the Galactic plane. The latitude profile of synchrotron 
intensity cannot be explained by a single emission component. 
The constancy of synchrotron intensity for b > 30° implies that 
the emission is originated from the Galactic halo. If it is com- 
ing from a thin disk, the intensity should scale as smb. If there 
is only the halo component, we will not have such a sharp in- 
crease of intensity as we approach to b = 0°. Therefore, we 
need to take into account additional thin disk component for 
the Galactic plane (i.e. b ^ 0°). 

Judging from visual inspection, the probability density func- 
tion (PDF) of synchrotron intensity for b > 30° can be ex- 
plained by a combination of two components. One component 
has a peak at T ^ 20A' and the other has a peak at T ^ 357*^. 
The latter may correspond to the emission from the North Po- 
lar Spur The former seems to have roughly a symmetric shape 
(Fig. O. However, the PDF of intensity times sinb shows a 
more complicated shape (dotted line). 

In this section, we study synchrotron emission from the 
Galactic halo (i.e. b > 30°) and the Galactic disk (i.e. \b\ < 2°) 
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separately. Our main goal is to see if statistics of synchrotron 
emission from the halo is consistent with turbulence models. 
When it comes to synchrotron emission from the Galactic disk, 
it is not easy to separate diffuse emission and emission from 
discrete sources. Therefore, we do not try to study turbulence 
in the Galactic disk. Instead, we will try to estimate which kind 
of emission is dominant in the Galactic disk. 

4.2. Structure of synchrotron halo 

How are the observed angular spectrum of synchrotron emis- 
sion and the spectrum of 3-dimensional MHD turbulence re- 
lated? In order to understand the relation, we need to under- 
stand the structure of the Galactic halo. There exist several 
models for the diffuse Galactic radio emission. Beuermann, 
Kanbach, & Berkhuijen (1985) showed that a two-component 
model, a thin disk embedded in a thick disk, can explain ob- 
served synchrotron latitude profile. They claimed that the 
equivalent width of the thick disk is about several kiloparsecs 
and thin disk has approximately the same equivalent width as 
the gas disk. They assumed that, in the direction perpendicu- 
lar to the Galactic plane, the emissivity e of each component 
follows 

e(z) = e{Q)sechiz/zo), (13) 

where z is the distance from the Galactic plane and zo is a con- 
stant. Recently, several Galactic synchrotron emission models 
have been proposed in an effort to separate Galactic compo- 
nents from WMAP polarization data (see, for example. Page 
et al. 2007; Sun et al. 2008; Miville-Deschenes et al. 2008; 
WaeUcens et al. 2008). All the models mentioned above assume 
the existence of a thick disk component with scale height equal 
to or greater than Ikpc. Sun et al. (2008) considered an addi- 
tional local spherical component motivated by the local excess 
of the synchrotron emission that might be related to the "local 
bubble". 

The detailed modeling of the Galactic synchrotron emission 
is beyond the scope of our paper We will simply assume that 
there is a thick component with a scale height of ~ Ikpc. We 
will also assume that there could be an additional local spherical 
component. Then, what will be the relation between spectrum 
of 3-dimensional turbulence and the observed angular spectrum 
of synchrotron emission? 

4.3 . Spectrum and Structure function of the 408MHz Haslam 

map 

In left panel of Fig.|4]we plot angular power spectrum of syn- 
chrotron emission at 408Mhz. The upper curve (solid curve) is 
the angular spectrum of the all-sky data. The lower curve (dot- 
ted curve) is that of high galactic latitude. We obtain the latter 
as follows. We first set the synchrotron intensity to zero for pix- 
els with \b\ < 30° and then calculate the angular spectrum using 
anafast packageQ Due to the partial sky coverage, the spectrum 
(lower curve) show an oscillatory behavior. The slope of the 
lower curve is very close to that of the straight line, which has a 
slope of -3. This result is consistent with earlier results for the 
408-MHz Haslam map (Tegmark & Efstathiou 1996; Bouchet, 
Gispert, & Puget 1996). Recently La Porta et al. (2008) per- 
formed a comprehensive angular power spectrum analysis of 

* HEALPix Homepage: |http://healpix.jpl. nasa.gQv/j 
' One should be careful if there is a white noise. In the presence of white noise, 
r(ei) — r(e2)P > + < |<5i -<52p >, where 5i and Si represent noise (Chepumov, 
of the slope. However, if lim9^o^2 is small enough, we can ignore the noise. Th 
is sufficiently smaller than values D2 shown in Fig.|4] 



all-sky total intensity maps at 408MHz and 1420MHz. They 
found that the slope is close to -3 for high galactic latitude re- 
gions. Other results also show slopes close to -3. For exam- 
ple, using Rhodes/HartRAO data at 2326 MHz (Jonas, Baart, 
& Nicolson 1998), Giardino et al. (2001b) obtained a - 2.92 
for high galactic latitude regions with \b\ > 20°. Giardino et 
al. (2001a) obtained a ^ 3.15 for high galactic latitude regions 
with \b\ > 20° from the Reich & Reich (1986) survey at 1420 
MHz. Bouchet & Gispert (1999) also obtained a ^ /"^ spectrum 
from the 1420 MHz map. 

Since we obtained angular spectrum from data with incom- 
plete sky coverage, the slope of the spectrum may be contami- 
nated by the edge effect. To avoid the effect, we may calculate 
angular correlation function K(6) first and then obtain C/ using 
Eq. ( fT9] l. But the angular spectrum Q obtained in this way is, 
in general, noisy. Therefore it is not easy to accurately measure 
the slope of the spectrum. When we want to obtain only the 
slope of the angular spectrum on small angular scales, hence 
the slope of the 3D spatial turbulence spectrum, we can use the 
second-order angular structure function: 

D2(0)=<|r(ei)-r(e2)P>. (14) 

As we discussed in ^3.2[ the slope of the second-order structure 
function will tell us about the slope of the angular spectrumQ. 

In the middle panel of Fig. |4] we show the second-order struc- 
ture function for the Galactic halo (i.e. \b\ > 30°). The slope 
of the second-order structure function lies between those of 
two straight lines. The steeper line has a slope of 4/3 and 
the other one has a slope of 1. The actual measured slope is 
^ 1.2. This result implies that the 3D turbulence spectrum is 
Ejoik) oc k~^-^, which is shallower than the Kolmogorov spec- 
trum £'3D(/t) oc r"/\ 

4.4. Effects of inhomogeneity and model calculations 

Then why is the slope shallower than that of Kolmogorov? 
Chepurnov (1999) gave a discussion about the effects of density 
stratification on the slope. He used a Gaussian disk model and 
semi-analytically showed that the slope of the angular spectrum 
can be shallower than that of Kolmogorov. In this subsection, 
we present further discussions. 

In the model discussed in |3] we assumed turbulence is 
homogeneous. However, this is certainly an unrealistic as- 
sumption for the Galactic halo. Synchrotron emission models 
(see previous subsection) assume either exponential (e""/"°) or 
square of hyperbolic secant {sech\z/zo\) law for synchrotron 
emissivity, where z is the distance from the Galactic plane. For 
simplicity, we assume the observer is at the center of a spherical 
halo. That is, the geometry is not plane-parallel, but spherical. 
In what follows, we use r, instead of z, to denote the distance to 
a point. 

To illustrate the effects of this inhomogeneity, we test 3 mod- 
els: 

1 . Homogeneous halo: Turbulence in halo, thus emissiv- 
ity, is homogeneous. Turbulence has a sharp boundary 
at dmax = ikpc. The outer scale of turbulence is lOOpc. 
Basically, this model is the same as the one we consid- 
ered in |3] 

second-order structure function will be Diid) =< |7'(ei) + 5i — r(e2) — <52p >=< 
private communication). Therefore, the noise can interfere accurate measurement 
is is true in our case. The Haslam 408MHz map shows D2(0.015°) ~ 0.05, which 
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Fig. 2. — Latitude distribution of average intensity. Tlie solid line shows the usual average intensity. The dashed line and the dotted line depict average taken near 
the Galactic center (—45° < Icai < 45°) and the Galactic anti-center (135° < Icai < 225°), respectively. Here Iq^i is the galactic longitude. 

Fig. 3. — PDF for \b\ > 30°. PDF of synchrotron intensity (solid curve) may be composed of two components: One has a peak at T 20K and the other at 
T ^ 35K. However, that of intensity time sin b (dotted curve) does not seem to be physically meaningful. 
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Fig. 4. — Haslam 408MHz map. Left: Angular spectra of full-sky and partial-sky (\b\ > 30°) maps. Spectrum of partial-sky map (thin curve) suffers from 
edge effect because pixels with \h\ < 30° are set to zero. Middle: Second-order angular structure function shows a slope of ^ 1.2, which is shallower than that of 
Kolmogorov turbulence (slope = 5/3). Right: Structure function as a function of galactic latitude. From bottom to top, the second-order angular structure functions 
are obtained for thin stripes {\Ah\ < 2°) along galactic latitudes of 30°, 10°, and 0°. 
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Fig. 5. — Halo models with stratification. Upper plot: Exponentially stratified halo. We take zo(= fo) = Ikpc and L (=eddy size) = lOOpc. Lower plot: Two- 
component halo. We take zi(= r[) = Ikpc, Z2(= ^l) = lOOpc, and L (=eddy size) = lOOpc. 



2. Exponentially stratified halo: Emissivity shows an ex- 
ponential decrease, e(r) oc e'''!''". We assume ro = Ikpc 
and the halo truncates at r = %kpc. The outer scale of 
turbulence is lOOpc. See Fig. |5] 

3. Two-component halo: Emissivity decrease as e(r) oc 
eie"''''''' +e2e~''/''-, where t2 = lOei, ri = Ikpc, and r2 = 



lOOpc. The halo truncates at r = 8kpc. The outer scale 
of turbulence is lOOpc. The second component mim- 
ics local enhancement of synchrotron emissivity. See 
Fig. 13 



We numerically calculate the angular correlation function 
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Fig. 6. — Model calculations. Three toy models for emissivity profiles in the Galactic halo are considered: homogeneous (solid line) halo, exponentially stratified 
halo (dotted line), and two-component exponential halo (dashed line). Left: Angular correlation functions. When angular separation is large, angular correlation 
functions follow the universal relation: (tt - S)/ sin S 1 /S. Middle: Second-order structure functions. When angular separation is small, the slope of the homo- 
geneous turbulence (solid line) is compatible with 5/3. But those of stratified halo models are shallower The two-component model (dashed Une) shows a slope 
compatible with -3, if we measure average slope between Q ~ 0.5° and ~ 10°. Right: Angular spectra. In general, the spectra are compatible with the 3D spatial 
turbulence spectrum of /""/3 for large /. The stratified halo models show shallower slopes, if we measure average slope between / ~ 10 and ~ 200 for instance. 
The homogeneous turbulence model also gives a shallower slope for these values of /. But, its spectrum shows a easily noticeable break near / ~ 50. 



K(9) and the second-order structure D2(0) from 

K(e) = J dn J dr2lC(\ri-r2\)e{n)e(r2), (15) 

D2{e)ccT-K(e), (16) 

where |ri-r2| = -I- - 2rir2Cos(?, e(r) is the synchrotron 
emissivity, T = limg^[)K(9), and we use the spatial correlation 
function /C(r) obtained from the relation: 

/C(r)oc / 4Trk^EiD(k)— — dk (17) 
Jo kr 

where the spatial spectrum of emissivity Eid has the form: 

/ constant if ^ < .le^ 
(k/ko^U^ ifk>ko, ^^^^ 

which is the same as Kolmogorov spectrum for k > ko ('^ 1 /L). 
The reason we use a constant spectrum for k < ko is explained in 
Appendix B (see also Chepurnov 1999). We obtain the angular 
spectrum from the relation: 

Q (X J P,(cose)K(cos0)dicos0), (19) 

where Pi is the Legendre polynomial. 

In Fig.|6] we plot the calculation results. The angular corre- 
lation function K(d) does not change much when 6 is small, and 
follows ^ (tt-O)/ sin9 ^ 1/6 when 9 is large. The critical an- 
gle is a few degrees for homogeneous model (thick solid curve) 
and single-component exponential model (dotted curve). As we 
discussed earlier, the critical angle for homogeneous turbulence 
is ^ (L/d„uuY"'' ^ 6°, where d,„ax (= Ikpc in our model) is the 
distance to the farthest eddy. In Fig.|6](left panel) we clearly see 
that the slope of K(6) changes near 0^6°. The second-order 
structure function D2{9) also shows a change of slope near the 
same critical angle {0 ^ 6°). In single-component exponential 
model (dotted curve), the value of d,„ax is not important. In- 
stead, the scale height zo is a more important quantity, which is 
\kpc in our model. In left and middle panels of Fig. |6] we ob- 
serve that the single-component exponential model also show a 
change of slope near ^ ?i few degrees. Therefore, we can in- 
terpret that the critical angle for stratified turbulence is ^ L/za, 
instead of ^ L/d„,ax 

Then, can we answer our earlier question of why the ob- 
served slope is shallower than that of Kolmogorov? Let us take 
a look at the right panel of Fig.|6l All 3 models show that the 



slope of c/ is almost Kolmogorov one for / > Icr ^ ndmax/L ^ 
30. However, if we measure average slope of c/ between / = 10 
and 200, we obtain slopes shallower than Kolmogorov. The 
single-component model and the homogeneous model give sim- 
ilar average slopes, which is a bit shallower than -3. However, 
the homogeneous model gives a more abrupt change of slope 
near I ^ 30. (In fact, right panel of Fig.|6]shows the break hap- 
pens near / ^ 50.) The average slope of the two-component 
model gives more or less gradual change of the slope and the 
slope is very close to -3 for a broad range of multipoles /. It is 
difficult to tell which model is better because the models are 
highly simplified. But, the two-component model looks the 
most promising. 

Note that, compared with the two-component model, the 
single-component model shows a more or less sudden change 
of slope near I ^ Icr ^ t^zo /L ^ 30. Therefore, if the single- 
component model is correct, the scale height zo cannot be much 
larger than ^10 times the outer scale of turbulence L. If zq 
is much larger than ^ lOL, Icr becomes smaller and we will 
have almost Kolmogorov slope for Z > 10. We also note that 
it is possible that 3D spatial turbulence spectrum itself can be 
shallower than the Kolmogorov one. That is, it is possible that 
spectrum of B(r), hence that of B^{r), can be shallower than 
the Kolmogorov one. For example, recent studies show that 
strong MHD turbulence can have a k~^-^ spectrum, rather than 
r"/3 (Boldyrev 2006; Beresnyak & Lazaiian 2006). If this is 
the case, the observed angular spectrum can be shallower than 
Kolmogorov for I > Icr- 

4.5. Synchrotron emission from Galactic disk 

In right panel of Fig. |4] we show how the second-order struc- 
ture function changes with galactic latitude. The lower curve is 
the second-order angular structure function obtained from pix- 
els in the range of 28° <b< 32°. The middle and upper curves 
are the second-order angular structure functions obtained from 
pixels in the range of 8° < < 12° and -2° < b < 2°, re- 
spectively. The middle and upper curves clearly show break 
of slopes near 6* ^ 3° and ^ 1 .5°, respectively. When the angu- 
lar separation is larger than the angle of the break, the structure 
function becomes almost flat. As we move towards the Galactic 
plane, the sudden changes of slopes happen at smaller angles. 

What causes this break of slope? There are at least two possi- 
bilities. First, a geometric effect can cause it. As we discussed 
in ^ change of slope occurs near 9c ^L/ dmax- As we move to- 
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wards the Galactic plane, the distance to the farthest eddy, dmax, 
will increase. As a result, the critical angle 9c ^ L/dmax will de- 
crease. Therefore, we will have smaller 9c towards the Galactic 
plane. This may be what we observe in the right panel of Fig. |4] 
Second, discrete synchrotron sources can cause flattening of the 
structure function on angular scales larger than their sizes. Al- 
though the map we use was reprocessed to remove strong point 
sources, there might be unremoved discrete sources. When fil- 
amentary discrete sources dominate synchrotron emission, the 
second-order structure function will be flat on scales larger than 
the typical width of the sources. In reality, both effects may 
work together At this moment, it is not easy to determine 
which effect is more important. 

4.6. On the polarized synchrotron emission 

Roughly speaking, the shape of the angular spectrum of po- 
larized synchrotron emission will be similar to that of the total 
intensity at mm wavelengths. However, at longer wavelengths, 
Faraday rotation and depolarization effects may cause flatten- 
ing of the angular spectrum, which has been actually reported 
(see de Oliveira-Costa et al. 2003 and references therein). On 
the other hand. La Porta et al. (2006) analyzed the new DRAO 
1 .4GHz polarization survey and obtained angular power spectra 
with power-law slopes in the range [-3.0,-2.5]. More observa- 
tions on the polarized synchrotron foreground emission can be 
found in Ponthieu et al. (2005), Giardino et al. (2002), Tucci et 
al. (2002), Baccigalupi et al. (2001). In this paper, we do not 
discuss further about polarized synchrotron emission. Read- 
ers may refer to recent models about the polarized synchrotron 
emission (Page et al. 2007; Sun et al. 2008; Miville-Deschenes 
et al. 2008; Waelkens et al. 2008). 

5. PROPERTIES OF THE DUST EMISSION MAP 

Thermal dust emission is also an important ingredi- 
ent of CMB foregrounds. In this section, we analyze 
a model dust emission map created by Finkbeiner et al. 
(1999), which is available at the NASA LAMBDA website 
(|Fttp://lambda. gsfc.nasa.gov/). 

In Fig. |7] we present statistical properties of the map. The 
map shows rough constancy of emission for high galactic lati- 
tude region when multiplied by sin/? (left panel of Fig. |7]). The 
sin/? factor also appears in the PDF: the sin/? factor makes the 
PDF more symmetric (middle panel of Fig. [T]). Therefore it is 
natural to conclude that a thin disk component dominates the 
dust map. 

Angular spectrum of the dust map shows a slope much flatter 
than that of the Haslam map (right panel of Fig. 13 . As in the 
Haslam map, the upper curve (solid curve) is the angular spec- 
trum of the all-sky data. The lower curve (dotted curve) is that 
of high galactic latitude. The slope of the dotted curve is very 
close to -2.5. However, due to the edge effect, it is not clear 
whether or not the slope is the true one. Since it is not affected 
by the edge effect, the second-order structure function can re- 
veal the true scaling relation on small angular scales. Note that 
the spectrum implies that the second-order structure func- 
tion has a slope very close to 0.5. Indeed, the actual second- 
order structure function of the dust map (see ® shows a slope 
of ^ 0.6, which corresponds to angular spectrum of ^ 

The slope of the angular power spectrum of the model dust 
emission map is very similar to that of the original FIR data. 
Schlegel et al. (1998) found a slope of -2.5 for the original FIR 
data. On the other hand, other researchers found slopes close to 



-3 from other observations (see Tegmark et al. 2000 and refer- 
ences therein; see also Masi et al. 2001). 

6. POLARIZED EMISSION FROM DUST 

Polarized radiation from dust is an important component 
of Galactic foreground that strongly interferes with intended 
CMB polarization measurements (see Lazarian & Prunet 2001). 
Then, what will be the angular spectrum of the polarized radia- 
tion from the foreground dust? One of the possible ways to esti- 
mate the polarized radiation from dust at the microwave range is 
to measure star-light polarization and use the standard formulae 
(see, for example, Hildebrand et al. 1999) relating polarization 
at different wavelengths. In this section, we describe how we 
can obtain a map of polarized emission from thermal dust and 
we discuss the angular spectrum of the polarized emission from 
thermal dust in the high galactic latitude (say, \b\ > 20°). 

6. 1 . Getting a map of polarized emission from thermal dust 

In principle we can construct a polarized dust emission map 
at mm wavelengths {Ipoi.mmil,b)) from a dust total emission map 
(Ii„m(l,b)) and a degree-of-polarization map {Pem,mm(l,b)) at mm 
wavelengths: 

^ pol jiim(J ^ b^ ~ Pem.mmiJ -jb^ ^mm(.^ :b\ (20) 

where (/,/>) denotes the galactic coordinate. However, neither 
Ipoi.mm{l,b) nor I,„„,{l,b) is directly available. Therefore, we 
need indirect methods to get Ipoi,mmil,b) and Imm(l,b). 

Obtaining a dust total emission map (/„„„(/,/>)) is relatively 
easy because dust total emission maps at FIR wavelengths are 
already available from the IRAS and COBE/DIRBE observa- 
tions. Using the relation 

I,n,„{l,b) = /ioo^,„(/,/')(1™«/100m'«)"'', (21) 
where 1 ^ /? < 2, one can easily obtain an emission map at 
mm wavelengths (/„„„) from the maps at 100 jim or 240 ^m. 
However, more sophisticated model dust emission maps at 
mm wavelengths akeady exist. For example, Finkbeiner et al. 
(1999) presented predicted full-sky maps of microwave emis- 
sion from the diffuse interstellar dust using FIR emissions maps 
generated by Schlegel et al. (1998). In fact, the model dust 
emission map we analyzed in the previous section (© is one 
of the maps presented in Finkbeiner et al. (1999). Therefore, 
we can assume that the thermal dust emission map (/„,„,(/,/')) is 
already available. 

Then how can we obtain a degree-of-polarization map at 
mm wavelengths {Pem.mm{l,b))l We can use measurements of 
starlight polarization at optical wavelengths to get Pem,mm(l,b). 
The basic idea is that the degree of polarization by emission at 
mm {Pem,mm) IS related to that at optical wavelengths {P em, optical), 
which in turn is related to the degree of polarization by absorp- 
tion at optical wavelengths (Pabs,opticai)'- 

Pahs, optical ^ Pem,optical ^ Pem.mm- (22) 

We describe the relations in detail below. 

When the optical depth is small, we have the following rela- 
tion (see, for example, Hildebrand et al. 2000): 

Pem.oplical ~ Pabs.oplical / 7 (23) 

where Pem.opiicai is the degree of polarization by emission and 
T is the optical depth (at optical wavelengths). We obtain po- 
larization by emission at mm wavelengths {Pem,min) using the 
relation 



em,mm ~ Pem,optical 



C -C 



Cmax Cm 



I 



Cm 



- Cm 



Cmax ~l~ Cm 



optical 



(24) 
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Fig. 7. — Average dust emission intensity times sini. Left: Latitude profile of averaged dust emission intensity time smb. Average value is roughly constant 
in high latitude region. The solid line represents the average taken over 360° . The dashed line and the dotted line depict average taken near the Galactic center 
(-45° < /go/ < 45°) and the Galactic anti-center (135° < la„i < 225°), respectively. Middle: PDF for \b\ > 30°. PDF of dust intensity times sinfo (dotted curve) 
shows a rough symmetry, but that of dust intensity (solid curve) does not. Right: Angular spectra of full-sky and partial-sky (\h\ > 30°) maps. Spectra are for dust 
map times sin fo. 
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Fig. 8. — The ratio of P(.m.mm//'(.ni.opi,ta;- Left- The polarization ratio vs. the axis ratio (132/01) of oblate spheroidal grains. The polarization ratio shows only a weak 
dependence on the axis ratio. We assume that the grain size is 0. l^m, magnetic field is perpendicular to the line-of-sight, \i,ptical = 0.5/.im, and A„„„ = 1000/xm. 
Right: The polarization ratio vs. the angle between magnetic field and the plane of the sky. We assume that the grain size is OA^m, the grain axis ratio is 1.5, 
Xoptical = 0.5/.im, and A„„„ = 1000/im. 



where C's are cross sections (of grains as projected on the sky) 
that depend on the geometrical shape (see, for example, the dis- 
cussion in Hildebrand et al. 1999; see also Draine & Lee 1984) 
and dielectric function e = e\ + i€2 (see Draine 1985) of grains. 

For mm wavelengths, it is easy to calculate the ratio in 
Eq. (l24b because the wavelength A is much greater than the 
grain size a (i.e. A ^ 2™). In this case, if grains are oblate 
spheroids with a\ < 02 = 03 and short axes (a\) of grains are 
perfectly aligned in the plane of the sky, we have 



A (L,.[ei(A)-l]+l) +[L,e2(A)]' 

where L values are defined by 

Li = [(l+/)//2][l-(l//)arctan/], 
L2 = L3 = (l-ii)/2, 

/ = (fl2/fli)'-l (26) 
(see, for example, Hildebrand et al. 1999). 



However, for optical wavelengths, the condition A :i> lira is 
not always valid and, therefore, the expression in Eq. dZSl l re- 
turns only approximate values. For accurate evaluation of the 
cross sections, one should use numerical methods. Fortunately, 
several numerical codes are publicly available for such calcu- 
lations (for example, DDSCAT package by Draine & Flatau 
1994, 2008; ampld.lp.f by Mishchenko 2000). We use am- 
pld.lp.f to calculate the ratio in Eq. ( l24b . We assume that 
the grains are oblate spheroids, grain size is 0.1//m, Xoptical = 
0.5/im, and A,,,,,, = lOOO/im. Left panel of Fig.[8]shows that the 
ratio is around 1 .5 when magnetic field is perpendicular to the 
line-of-sight. It also shows that the ratio of Pem,mm/Pem.opiicai is 
almost independent of the grain axis ratio. 

In this subsection, we described a simple way to obtain a po- 
larized map at mm wavelengths. However, actual implementa- 
tion of the method can be more complicated due to the follow- 
ing reasons. First, we used an assumption that the grains that 
produce optical absorption produce also microwave emission. 
But, this is not true in general (see Whittet et al. 2008). Second, 
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the expressions in Eqs. ( l24b and dZSl ) are vahd when magnetic 
field direction is fixed and perpendicular to the line-of-sight and 
all grains are perfectly aligned with the magnetic field. If this 
is not the case, Eq. jTM will become Pem.mm oc Pem.opiicai with 
the constant of proportionality that depends on magnetic field 
structure and the degree of grain alignment. The effect of partial 
aUgnment is expected to be less important Then, how serious 
is the effect of non-perpendicular magnetic field? We perform 
a numerical calculation using ampld.lp.f to evaluate the effect. 
We assume that the grains are oblate spheroids, grain size is 
0.1/im, Xoptkai = 0.5/im, and A„,„, = 1000/im. Right panel of 
Fig. [8] shows that the polarization ratio drops from ^ 1.5 to 
^1.1 when the angle (between magnetic field and the plane of 
the sky) changes from 90° to ~ 5°. Therefore, the effect is not 
very strong and can be potentially corrected fofl 

6.2. Angular spectrum of polarized emission from thermal 

dust 



After we have constructed a map of the polarized emission 
from thermal dust, we can obtain the angular spectrum. How- 
ever, if we are interested in the shape of the angular spectrum, 
we do not need to construct the polarized thermal dust emission 
map. We can get the shape of the the angular spectrum directly 
from the starlight polarization map Pabs.opticai{l,b). 

Eq. (HOll tells us that /,;„/,„„„ is given by Pem.mm times /,„„,. 
From Eqs. ( |23] i and ( |24] |. we have 

ipoLmm ~ Pem^mm hnm ^ Pem^optical hnm 
~ iPabs,optical Z"^) ^mm 

^ Pabs .optica! • (27) 

Here we use the fact r oc /„„„. Note that the constant of propor- 
tionality does not affects the shape of the angular spectrum if 
grain properties do not vary much in halo. Therefore, as to the 
power spectrum c/ of Ipoi.mm, we can use that of Pabs.opiicai- 

Ci of Ipol .mm ^ Cl Pabs, optical- 

(28) 

Once we know the angular spectrum of Pahs,opticai, we can esti- 
mate the angular spectrum of Ipoi.mm- 

Then, what is the measured angular spectrum of 
Pabs.opiicaii^,^)'^ Fosalba ct al. (2002) obtained C/ ~ /"''^ for 
starlight polarization. The stars used for the calculation are at 
different distances from the observer and most of the stars are 
nearby stars. The sampled stars are mostly in the Galactic disk. 
CL02 reproduced the observed angular spectrum numerically 
and showed that the slope becomes shallower when only stars 
with a large fixed distance are used for the calculation. There- 
fore, it is clear that distance, or dust column density, to the stars 
is an important factor that determines the slope. We expect that, 
if we consider only the nearby stars with a fixed distance, the 
slope will be steeper This means that, if we consider stars in 
the Galactic halo, the slope will be steeper 

The method described above requires measurements of po- 
larization from many distant stars in the Galactic halo. Unfortu- 
nately, the number of stars outside the Galactic disk that can be 

^ Theories (see Lazarian 2007 for a review) predict tliat grains starting with a particular size get aligned. If the grain size distribution varies from one place to another, 
this fraction of aligned grains will also vary. However, for the diffuse ISM, the grain distribution should not vary much (Weingartner & Draine 2001). Therefore, the 
effect of partial alignment will be less significant. 

' We can make use of Right panel of Fig.[8]reversely. In the future, when we can accurately measure polarized emission from thermal dust in FIR or mm wavelengths, 
we can obtain the values of [(Cmax-Cmm)/(Cmax + C,„i„)]i„m- This result combined with the values of [(C,„ax-C„i„)/{C,„ax + Cmin)]optical ™ optical wavelengths can 
be used to find average angle between magnetic field and the plane of the sky. That is, when we know the ratio [■■■]mm/[---]opiicah we can use Right panel of Fig.[8]to 
find the angle between magnetic field and the plane of the sky. 

Consider a 3-dimensional magnetic field with a 3D spectrum Ej[)(kx,ky,k^). The spectrum of the magnetic field on a two-dimensional sub-plane (e.g. z = plane) 
is E^^piane(kx,ky) OC dk^ E-i[)(kx,ky,kz), which we use in our calculations to generate two sets of magnetic field on a two-dimensional plane. 



used for this purpose are no more than a few thousands (Heiles 
2000; see also discussions in Page et al. 2007; Dunkley et al. 
2008). When more observations are available in the future, ac- 
curate estimation of Ipoi,mm{l,b) (and C/ of Pabs,opticai) will be 
possible. We do not pursue this topic further in this paper. 

6.3. Model calculations for starlight polarization 

Then, what do we expect about C/ of Pabs.opiicai for the Galac- 
tic halo? To deal with this problem we use numerical simula- 
tions again. We first generate two sets of magnetic field on a 
two-dimensional plane (8192 x 8192 grid points), using Kol- 
mogorov three-dimensional spectrcQ We consider 3 models: 

1. Case 1, Nearby stars in a homogeneous turbulent 
medium: We generate three (i.e. x,y, and z) components 
of magnetic field on a two-dimensional plane (8192 x 
8192 grid points representing AQQpc x AQQpc), using 
the following Kolmogorov three-dimensional spectrum: 
E^oik) oc r'l/^ \fk> ko, where fc,, - 1/100 pc. (The 
outer scale of turbulence is lOOpc.) We assume the vol- 
ume density of dust is homogeneous. All stars are at a 
fixed distance of lOOpc from the observer. 

2. Case 2, Distant stars in a homogeneous turbulent 
medium: We generate three (i.e. x,y, and z) components 
of magnetic field on a two-dimensional plane (8192 x 
8192 grid points representing 4 kpc x 4 kpc), using 
the following Kolmogorov three-dimensional spectrum: 
Ej,D(k) oc r"/^ ifk> ko, where ko - 1 /lOO pc. Other 
setups are the same as those of Case 1, but the distance 
to the stars is 2kpc. 

3. Case 3, Stars in a stratified medium: We use the mag- 
netic field generated in Case 1 . The volume density of 
dust shows a sech^(z) decrease: p(r) = 4-po/[e\p(r/ro) + 
exp(— r/ro)]^. We assume spherical geometry and ro = 
lOQpc. The stars are at r = 2QQpc from the observer. 
The outer scale of turbulence is lOOpc. 

We assume that dust grains are oblate spheroids. In the pres- 
ence of a magnetic field, some grains (especially large grains) 
are aligned with the magnetic field (see Lazarian 2007 for a re- 
view). Therefore, cross sections parallel to and perpendicular to 
the magnetic field are different. We assume that parallel cross 
section is ^30% smaller than the perpendicular one. We use 
the following equations to follow changes of Stokes parameters 
along the path: 



ndI/ds = -5 + AaQ/I, (29) 

d{Q/I)/ds = Aa-A<j(Q/I)\ (30) 

d{U/I)/ds = -Aa(Q/I){U/I) (31) 

(see Martin 1974 for original equations; see also Dolginov, 

Gnedin, & Silantev 1996), where 6 = (cti + 172), Act = (cti -CT2), 
and 

2ai=a±, (32) 

2ct2 = ctj_-(ct^-CT||)cos7 (33) 
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(Lee & Draine 1985). Here tr^ and a\\ are the extinction coef- 
ficients and 7 is the angle between the magnetic field and the 
plane of the sky. After we get the final values of Stokes parame- 
ters, we calculate the degree of polarization (•\/ Q^ + U^ /I) and, 
then, the second-order angular structure function of the degree 
of polarization. 

We show the result in Fig. |9] When the stars are nearby 
(Case 1), the spectrum is consistent with the Kolmogorov spec- 
trum for small 6. The result for the stratified medium (Case 
3) also shows a spectrum compatible with the Kolmogorov one 
for small 6. When stars are far away (Case 2), the qualitative 
behavior is similar. However, if we measure average slope be- 
tween 9 = 0.2° and 20°, the result is different; the slope for 
Case 2 is substantially shallower Note that 9 = 0.2° and 20° 
correspond to / = 1000 and 10, respectively. This means that, 
when we have either distant stars or mixture of distant and 
nearby stars, the angular spectrum will be shallower than the 
Kolmogorov one. Therefore, it is not surprising that Fosalba 
et al. (2002) obtained a shallow spectrum of ~ Q oc Z"'^ for 
stars mostly in the Galactic disk. Flattening of spectrum (i.e. 
Ci oc with a « 1.3 1.4) for polarized FIR dust thermal 
emission is also observed in Prunet et al. (1998; see also Prunet 
& Lazarian 1999). 

6.4. Comparison with the CMB polarization 

On the right panel of Fig. |9] we plot angular power spectrum 
of star light polarization. As we mentioned earlier, the angular 
spectrum of the degree of starlight polarization should be simi- 
lar to that of polarized thermal dust emission (see Eqs. 1271 and 

Esl). 

To obtain angular spectra, we use a Gauss-Legendre quadra- 
ture integration method as described in Szapudi et al. (2001). 
To be specific, we first generate magnetic fields from the 3 mod- 
els we considered in the previous subsection. Then, we calcu- 
late angular correlation functions, K{cos9). Finally, we obtain 
the angular spectra using Eq. ( fT9] l. Since C/ obtained in this 
way is very noisy, we plot C/ averaged over the multipole range 
(//1. 09, 1.09/). We do not show Q for / > 1000 because it is 
too noisy even with the averaging process. We normalize the 
spectra using the condition 1(1+ l)C//27r = 3(pK^). This 
normalization is based on the values given in Page et al. (2007; 
their Eq. [25]). We assume that the observed band is W-band 
(f = 94GHz). 

The plot shows that the slopes for nearby stars (thick solid 
line) and stars in the stratified medium are shallower than that of 
the Kolmogorov spectrum for / < 1000. This result is consistent 
with that obtained with the angular structure function. Note that 
the case of distant stars has much flatter spectrum for I < 1000. 
We believe that our toy model for the stratified medium (Case 
3) better represents the actual situation for polarized emission 
from thermal dust in the Galactic halo. Therefore, we expect 
that the polarized thermal emission from thermal dust in high- 
latitude Galactic halo has a spectrum slightly shallower than the 
Kolmogorov spectrum for / < 1000. Our calculation do not tell 
us about the slopes for / > 1000. However, judging from the 
behavior of the angular structure function for 6* < 0.1°, we ex- 
pect that Ci cx r'^/^ for I > 1000 (see the straight dashed Hne 
for / > 1000 on the right panel of Fig. |9l). 

The Figure shows that the EE spectrum dominates polar- 
ized thermal dust emission from high-latitude Galactic halo for 

' ' Even in the case of Planck checking the consistency of the foreground with the 
test, however. 



I > 100. The EE spectrum is expected to be sub-dominant when 
/ > 5000. 

7. UNDERLYING SPECTRUM: ADDITIONAL INFORMATION 
FOR REMOVING FOREGROUNDS 

Removal of Galactic foregrounds has always been a big con- 
cern for CMB studies. The challenge is only going to in- 
crease substantially now, when CMB polarization studies are 
attempted. 

The knowledge that the foreground are not an arbitrary noise, 
but have well defined statistical properties in terms of their spa- 
tial power spectra is an important additional information that 
can be utilized to evaluate and eventually eliminate the fore- 
ground contribution. 

Utilizing the information about underlying turbulence power 
spectrum is not straightforward, however. Our study shows that 
the observed power spectrum may depend on geometry of the 
emitting volume. Therefore, the detailed modeling of the fore- 
ground fluctuations should involve accounting for the geometry 
of the emitting volume. 

The latter point stresses the synergy of the Galactic fore- 
ground and CMB studies. Indeed, our fitting of the power spec- 
tra in Fig. |6] shows that on the basis of its variations we may 
distinguish between different models of the emitting turbulent 
volume. As soon as this achieved, one can predict, for instance 
the level of fluctuations that are expected from the foreground 
at the scales smaller than those studied. A simplification that 
is expected at higher resolutions that are currently available, is 
that at sufficiently small scales the statistics should get indepen- 
dent of the large-scale distributions of the emitting matter. 

Consider an example of utilization of this approach. As we 
know, Planck will sample the polarized foregrounds with the 
resolution of up to / ^ 2000 (69 ~ 5')- For Planck high fre- 
quency coverage should help to remove the foregrounds with 
high accurac}03. A higher resolution will be available, for in- 
stance, with balloon-borne missions, which in many cases will 
not have as many frequency channels as Planck does. How- 
ever, our present work suggests that the contribution of the 
foregrounds to the spectrum measured by these high-resolution 
missions can be evaluated on the basis of our knowledge of the 
underlying spatial power spectrum. 

Consider, for instance, Fig.|9l If Planck measures the spatial 
spectrum up to / = 2000, then for a higher resolution balloon 
mission one can evaluate the level of foreground contamination 
by extrapolating the expected foreground spectrum. 

While most of the paper is directly related to making use of 
the knowledge of the underlying spectra and/or two points cor- 
relations in order to deal with the foregrounds, the part dealing 
with higher-order statistics is not directly related to the fore- 
ground removal. Nevertheless, our analysis shows that high- 
order correlations can also help in understanding of the proper- 
ties of foregrounds. 

8. HIGH-ORDER STATISTICS 

High-order structure functions are used for the study of inter- 
mittency, which refers to the non-uniform distribution of struc- 
tures. The structure functions of order p for an observable / is 
defined by 

Sp(r)=<\I(x)-I(x + r)\P>, (34) 
lei of emission from a turbulent medium with a given geometry may be a useful 
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Fig. 9. — The second-order angular structure function and angular power spectrum for starlight polarization. Plots are for degree of polarization. Left: The 
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for the dashed line, which we believe is the most realistic model among the 3 toy models. We assume that the observed frequency is 94GHz. Note that 6 = 0.1° 
corresponds to / ~ 180°/6*° ~ 1800. We also show the polarized CMB 'EE' spectrum. 



where the angled brakets denote average over position x. For an 
observable defined in the plane of the sky, the angular structure 
function of order p is 

D^(0)=<|/(ei)-/(e2)|''>, (35) 

where Ci and 62 are unit vectors along the lines of sight and 9 is 
the angle between Ci and 62- 

Traditionally, researchers use high-order structure functions 
of velocity to probe dissipation structures of turbulence. In 
fully developed hydrodynamic turbulence, the (longitudinal) 
velocity structure functions Sp =< ([v(x + r)- v(x)] •f)'' >=< 
(5v^(r) > are expected to scale as r^f. One of the key issues 
in this field is the functional form of the scaling exponents (p. 
There are several models for (p. Roughly speaking, the dimen- 
sionality of the energy dissipation structures plays an important 
role. 

Assuming 1 -dimensional worm-like dissipation structures. 
She & Leveque (1994) proposed a scaling relation 

Cf = W9 + 2[l-(2/3)''/^] (36) 

for incompressible hydrodynamic turbulence. On the other 
hand, assuming 2-dimensional sheet-like dissipation structures, 
Miiller & Biskamp (2000) proposed the relation 

Cr = /7/9+l-(l/3)''/3 (37) 

for incompressible magneto-hydrodynamic turbulenc43 

Recently, high-order structure functions of molecular line in- 
tensities have been also employed (Padoan et al. 2003; Gustafs- 
son et al. 2006). In optically thin case, the molecular line inten- 
sities are proportional to the column density. Kowal, Lazarian, 
& Beresnyak (2007) studied scaling of density fluctuations in 

Boldyrev (2002) obtained the same scaling relation for highly supersonic turbulence. However, since it is unlikely that turbulence in the Galactic halo is highly 
supersonic, we refer the scaling relation to the "Muller-Biskamp" scaling. 

Consider a circular uniform cloud with radius A and intensity / centered at the origin. (For simplicity, let us consider the 2D Cartesian coordinate system.) Then 
we can show that D„(r) oc /"(27rr)A for r 2> A. When r ~ A, r-dependence will be weaker. We can show that structure functions for a filament also show a similar 
r' -dependence, because a filament can be viewed as a chain of circular clouds (or a chain of square-like clouds): D„(r) oc I"(2nr)LfA, where Lf is the length of 
the filament. Therefore, both filaments and cloud-like sources exhibit r'-scahng. The high-order structure functions for the dust emission map show ~ r' power-law 
scaling for 6^1° (left panel of Fig. I lit . Thus the dominant discrete sources are either thin filaments or point sources. We expect the typical size or width of the 
dominant structures is < 1 ° . 



MHD turbulence. Their numerical results show that behavior 
of the scaling exponents for column density depends on sonic 
Mach number of turbulence. 

The Haslam map shows a reasonable agreement with the 
Muller-Biskamp MHD model. In some sense this is natural 
because the synchrotron emission arises from MHD turbulence 
in the Galactic halo and roughly reflects the column density of 
cosmic-ray electrons. However, it is not clear why the column 
density shows a similar scaling as velocity. Note that Padoan 
et al. (2003) also obtained a similar result using ^^CO emission 
from Perseus and Taurus. However, the physical origin may be 
different: their result may reflect existence of highly supersonic 
turbulence in those clouds. 

Unlike the Haslam map, the dust map does not show agree- 
ment with existing models. The scaling exponents do not show 
strong dependence on the order p. Left and middle panels of 
Fig. [TT] show that the slope is around 1 for high-order struc- 
ture functions. This kind of behavior is expected when discrete 
structures dominate the map. However, it is not clear what 
kinds of discrete structures dominate. Either thin filamentary 
structures or point sources could explain the behavior of the 
structure function^ Since the Haslam map and the dust map 
sample different types of the ISM, it is not so surprising that 
they show different scaling behaviors. 

For the Galactic disk, high-order structure functions of both 
maps show nearly flat structure functions. The structure func- 
tions show a nearly flat behavior even for 6* ^ 1°, which means 
that the typical size or width of the dominant structures is ~ 1° 
or smaller. Again, this kind of behavior is expected when thin 
discrete filamentary structures or point-like sources dominate in 
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the Galactic disk. Note that, since the region of the sky we con- 
sidered is only a thin stripe (-2° < b < 2°) along the Galactic 
disk, the calculation domain becomes 1D-Iik43 and it is possi- 
ble that thin filaments and point sources of size < 1 ° can exhibit 
a slope flatter than r' . 

9. DISCUSSION AND SUMMARY 

In this paper, we have discussed the relation between 3D spa- 
tial spectrum and the observed angular spectrum (or the second- 
order angular structure function) of the CMB foreground emis- 
sions. We have focused on synchrotron total intensity and po- 
larized thermal dust emission. Our current study, as well as 
earlier studies (Chepurnov 1999; CL02), predicts that C; will 
reveal true turbulence spectrum on small angular scales. Then, 
on what scales will we see the turbulence spectrum? Our model 
calculations that take into account stratification effects imply 
that 

1. 9 < a few times 0.1° (or Z > a few times 100) for syn- 
chrotron emission (see Fig.|6]l, and 

2. 0<O.l°(or/>l 000) for polarized emission from ther- 
mal dust (see Fig.|9l). 

On larger angular scales, spectra are expected to be shallower. 

Then, how is the turbulence spectrum related to the angular 
spectrum (or the second-order angular structure function) on 
small scales? When the 3D spatial spectrum has a power-law 
spectrum, k~"\ the observed angular spectrum and the second- 
order angular structure function will be 

Ci (X r"' and (38) 
D2(0)cxr-2 if0«L/fl'™,„ (39) 

where L is the outer scale of turbulence and d,nax is the distance 
to the farthest eddy in case of homogeneous turbulence. We can 
interpret d„,a.i ~ zo if there is a stratification with scale height zo 
(see Fig. |6l). When the angular separation is large, we have a 
universal scaling for the angular correlation function: 

K(e)^^^l if 9 >L/d,„a, (40) 
smt^ 9 

In this limit, the structure function is roughly constant and the 
angular spectrum is roughly proportional to . 

In this paper, we have analyzed Haslam 408Mhz map and a 
model dust emission map and compared the results with model 
calculations. We have found that 



1 . The Haslam map for high galactic latitude (b > 30°) can 
be explained by MHD turbulence in the Galactic halo. 
The measured second-order angular structure function 
is proportional to 0'^, which corresponds to an angu- 
lar spectrum of T^^. The high-order statistics for high 
galactic latitude (b > 30°) is consistent with that of in- 
compressible magnetohydrodynamic turbulence. Our 
model calculations show that a two-component model 
(see 34.41 and Fig.|6]l can naturally explain the observed 
angular spectrum. The one-component model can also 
explain the observed slope. But, the slope of the spec- 
trum shows a more abrupt change near / ~ 30. 

2. The model dust emission map may not have anything to 
do with turbulence on large angular scales. That is, we 
do not find signatures of turbulence in the map. 

3. Both maps show flat high-order structure functions for 
the Galactic plane. This kind of behavior is expected 
when discrete structures dominate the map. 

We have described how we can obtain angular spectrum of 
polarized emission from thermal dust in high galactic latitude 
regions. Our model calculations show that starlight polariza- 
tion arising from dust in high galactic latitude regions will have 
a Kolmogorov spectrum, C/ oc /""/^, for / > 1000 and a shal- 
lower spectrum for / < 1000 (Fig|9]l. We expect that polarized 
emission from the same dust also has a similar angular spec- 
trum. That is, we expect that the angular spectrum of polarized 
emission from thermal dust is close to a Kolmogorov one for 
I > 1000. 
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APPENDIX 

A. THE SECOND-ORDER ANGULAR STRUCTURE FUNCTION IN THE SMALL ANGLE LIMIT 

The angular correlation K{9) is given by the integral 

K{e) = Jdh j dhJC{\h-h\), (Al) 

where JCif) is the 3 dimensional spatial correlation. Suppose that li is along x-axis, li = (/i,0), and I2 = Z2(cos6',sin6'). Then, the 
correlation in the limit of small 9 is given by 

' ~ (A2) 

(A3) 



j dh 


1 dh 


< 
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~ JdhJ dh IC {^{h-h)^ + hhe^) . (A4) 

Suppose that the spatial correlation follows a power law: /C(r) oc const -r^ for r < L, where L is the outer-scale of turbulence. For 
Kolmogorov turbulence, p = 2/3. Then the derivative of K{d) is given by 

^ oc -Jdh j dh [(h-hf + 1,126^]" '^~\2hh9) (A5) 

-oc Jdujdw [v? + {u^-vp-)e^/AY'^~\u^-v?)e/2, (A6) 

where u = l\+l2 and w = l\ — l2- \f p <\, the integration diverges as 9 goes to zeroH Therefore, when p < 1, it suffices to perform 
the integration in the vicinity of l\ = I2 or w = 0. Then we have 

^oc- Jdujdw [w^ + u^ey^Y^^-' (A7) 
Ki- J du(ue/2y-^u^9/2(x-9P, (A8) 

where we use dw/{w^+A^)" = A^~^" J^^j^ dO sec^~^" 9. Therefore, for small 6 we have 

K(e)(xCi-C20P^\ (A9) 
where Ci and C2 are constants. Comparing this equation with 

K(e) = C3- €402(6), (AlO) 

we get 

D2(0)(yi9''^K (All) 

Analytic expressions for the relation between the angular structure function {D2) and the spatial ID spectrum (E(k); in case of 
Kolmogorov, E(k) oc r^/^) can be found in the literature. For example, Lazarian (1995a; see also Lazarian & Shutenkov 1990) 
derived the following expression; 

/■^/* d 

E{k)^k dr]—{Q{r])Tj)MkRf]) + ICs, (A12) 
Jo drj 

where C can be regarded as the outer scale of turbulence, R is the size of the system, Q{rj) ^ 02(77)77, 77 = sin^, 7i(x) is the Bessel 
function of the first order, and IC5 is a small correction term. 

B. SPATIAL SPECTRUM OF EMISSIVITY 

The synchrotron emissivity is proportional to ~ n(e)B'' oc B^, where n{e) is the high-energy electron number density. Suppose that 
magnetic field is roughly a Gaussian random variable. This may not be exactly true, but should be a good approximation. When a 
Gaussian random variable B(r|3 follows a Kolmogorov spectrum 

EB.o^\Bikt^[l,,^^yn,. llll^ (Bl) 

we can show that the 3D spectrum of B^{r) follows Eq. [T8](see, for example, Chepurnov 1999). The correlation of B^{r) and 3D 
energy spectrum of B^{r) are related by 

/CB2(r)=<B^(x)B2(x + r)>,oc J Es2^^Di^)e"' ''d\ (B2) 

Es23d(^) = \BHkt « J ICB2{r)e-''' 'd\ (B3) 

where < ... >^ denotes an average over x. A Gaussian random variable satisfies 

< B^{x)B^(x + r) >=< B^(x) >< ^^(x + r) > +2 < B(x)B(x + r) >^ (B4) 

where the first term on the right is a constant. Therefore we can ignore the term in what follows. Fourier transform of both sides 
results in 

LHS=£B^3D(k), (B5) 

RHS = 2 [ < B(x)B(x + r) >^ g-^ 't/^r (B6) 



" When p= I, the spatial correlation becomes JC(r) oc C-r, where C is a constant. The corresponding 3D spectrum is E{k) oc k'^. When the slope of the turbulence 
spectrum is steeper than k^"*, the correlation function has the form K.(r) oc A2o-r' regardless of the turbulence slope. On the other hand, when the three-dimensional 
spectrum of turbulence is shallower than k^, we have IC{r) oc /Cq — r'""', where Ko ~ L'"^^ is a constant. Therefore, the condition of p < 1 is generally satisfied in 
turbulent medium. 

For simplicity, we assume i? is a scalar. 
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= 2 J d\ /CsCr) /CB(r)e"*-' (B7) 

= 2jd\j d^p j t/^q£B,3D(p)£B,3D(q)e'^''"''' '' (B8) 

= 2j d'pj d^q£B,3o(p)£B,3Z)(q)^(P + q-k) (B9) 

= 2j d^pEBMP)EB,3DQs.-P), (BIO) 
where 6{k) is the Dirak (5-function. Therefore we have 

£B^3Z)(k) « Eg^^iom « 2 y d\\EBM^)? ~ constant (B 1 1) 

for k <^ ko. 



